Phase Space Interpretation of Exponential Fermi Acceleration 
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Recently, the occurrence of exponential Fermi acceleration has been reported in a rectangular 
billiard with an oscillating bar inside [K. Shah, D. Turaev, and V. Rom-Kedar, Phys. Rev. E 
81, 056205 (2010)]. In the present work, we analyze the underlying physical mechanism and show 
that the phenomenon can be understood as a sequence of highly correlated motions, consisting of 
alternating phases of free propagation and motion along the invariant spanning curves of the well- 
known one-dimensional Fermi-Ulam model. The key mechanism for the occurrence of exponential 
Fermi acceleration can be captured in a random walk model in velocity space with step width 
proportional to the velocity itself. The model reproduces the occurrence of exponential Fermi 
acceleration and provides a good ab initio prediction of the value of the growth rate including 
its full parameter-dependency. Our analysis clearly points out the requirements for exponential 
Fermi acceleration, thereby opening the perspective of finding other systems exhibiting this unusual 
behaviour. 

PACS numbers: 05.45.a,05.45.Ac,05.45.Pq 
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I. INTRODUCTION 

Recently, the investigation of Fermi acceleration (FA) 
in two-dimensional (2D) time-dependent billiards has at- 
tracted a lot of attention [l|-[7J . Fermi acceleration is the 
unbounded energy gain of particles exposed to driving 
forces and was first proposed in 1949 by Enrico Fermi [8] 
to explain the high energies of cosmic radiation (for a re- 
view of FA see Ref. @). He suggested that charged par- 
ticles repeatedly interact with time-dependent magnetic 
fields (originating either from shock waves of supernovae 
or from magnetized interstellar clouds) in such a way that 
on average they gain energy. Nowadays, FA is investi- 
gated in a variety of systems bel ong i ng to different areas 
of physics, such as astrophysics |10l - [l2l ]. plasma physics 
|13l ll4J . atom optics [ljl [l6| and has even been used 
for the interpretation of experimental results in atomic 
physics J17j |. 

The one-dimensional (ID) prototype system allowing 
the investigation of FA is the so-called Fermi-Ulam model 
(FUM) [lj|, which consists of non-interacting particles 
moving between one fixed and one oscillating wall. The 
FUM and its variants have been the subject of extensive 
theoretical (see Ref. [18|, [l9| and references therein) and 
experimental [2fJ, [21| studies. In the FUM, the existence 
of FA depends exclusively on the driving-law of the os- 
cillating wall: As long as the driving law is sufficiently 
smooth, there is no unlimited energy growth due to the 
existence of invariant spanning curves [18]. In particular, 
this means that harmonic driving laws do not lead to FA 
in the FUM. 
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In 2D time-dependent billiards, already a smooth driv- 
ing law may lead to FA. For example, the existence of FA 
was shown for a harmonically oscillating stadium-like bil- 
liard [3j, |4| , in the driven eccentric annular billiard 5] , in 
an oval billiard [6[1 and in the time-dependent elliptical 
billiard p], 122, [231 . On the other hand, the breathing 
concentric annular Q and the circular billiard [2J| do 
not exhibit FA. The ensemble-averaged energy E(t) in 
all the 2D time-dependent billiards that do show FA, 
grows as a function of time according to a power law, 
i.e. E{t) ~ t d , with some exponent d. For such power 
laws, it is known that FA is not structurally stable in the 
sense that any finite amount of dissipation will destroy 
it [l|, 12514271 ] , independent on whether dissipation is in- 
troduced via inelastic collisions or via drag forces (e.g. 
Stokes' friction). 

Thus, a natural question to ask is whether there are 
certain time-dependent billiards that show a somewhat 
'faster' acceleration of the energy, in particular whether 
there is e.g. exponential acceleration possible. The first 
hint that such a fast acceleration process is possible is 
given in Ref. [28j , where the authors prove the existence 
of single orbits with exponential energy growth under cer- 
tain conditions. In Ref. 29] it is shown that the energy 
of a whole ensemble of particles grows exponentially in a 
rectangular billiard with an oscillating bar inside. This 
result is generalized in Ref. [2|, where the authors show 
by means of an analysis of the Anosov-Kasuga invariant 
that in special classes of billiard systems the ensemble 
averaged energy accelerates exponentially. 

While in Refs. @, [28|, [29[ a mathematical analysis of 
exponential FA is provided, the aim of the present work 
is to investigate the physical mechanism leading to expo- 
nential FA in the setup proposed in Ref. [29[ . This setup 
consists of an oscillating bar inside a rectangular billiard, 
where the bar is aligned parallel to the long side of the 
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FIG. 1. (Color online) Setup: Rectangular billiard of height 
H and length L with an oscillating bar of length I placed in 
the center, parallel to the a;-axis (here H = 2, L = 4, 1 = 2). 
The two areas above and below the bar are the interactions 
areas (with the moving bar) A and B, respectively. 



rectangle. This can be interpreted as particles moving 
alternately in a FUM and in a static rectangular billiard. 
Since neither the static rectangular billiard nor the FUM 
alone even show FA, we want to clarify from a physical 
point of view how the combination of the both leads to 
exponential FA, i.e. what are the microscopic processes 
that cause the astonishingly fast acceleration. To this 
end, we will show that in the high velocity regime, the 
temporal movement of the particles on invariant curves of 
the FUM can be modeled by a suitable random walk with 
step sizes being proportional to the velocity itself. This 
random walk model shows exponential acceleration. Fur- 
thermore, the corresponding parameters of the random 
walk can be extracted from the underlying FUM, even 
enabling an alternative (compared with the one given in 
Ref. [29]) prediction of the exponential acceleration rate 
without any free parameters. 

The work is structured as follows: In section [TT| we 
introduce the setup and show the results of our numeri- 
cal simulations. How the microscopic dynamics of single 
trajectories can be interpreted as piecewise motion along 
invariant spanning curves of an appropriate Fermi-Ulam 
model is shown in section ILUl With this picture in mind, 
we construct a random walk model in section IIVI that 
generalizes the considered model incorporating the ba- 
sic characteristics of the underlying physical mechanism. 
Finally, a short summary and outlook is given in section 

El 



II. SETUP AND RESULTS 

The investigated setup is shown in Fig. Q] It consists 
of a rectangular billiard of length L and height H with 
an oscillating bar inside. The rectangular billiard with- 
out the bar is integrable, since upon collisions with the 
billiard boundary, just the sign of the corresponding com- 
ponent of the velocity v is reversed, i.e. in particular \v x \ 



and \v y \ are preserved. Now, a bar of length I is placed in 
the middle of the billiard, parallel to the longer side (in- 
direction) of the billiard. Here, we assume an harmonic 
oscillation law: j/f, = acos(wi), with a and u) being the 
driving amplitude and frequency, respectively. Since the 
oscillating bar transfers momentum in y-direction only, 
\v x \ is preserved. As an initial ensemble, we take N = 10 4 
classical, non-interacting particles with a fixed velocity 
v x = 0.16; and v y is randomly chosen in the interval 
[0,40vJ (we use H = 2, L = 4, I = 2 for the simula- 
tions). We iterate these particles by numerically solving 
the corresponding discrete mapping, i.e. by calculating 
the successive collisions with the billiard boundary, which 
consists of the rectangle and the bar. The main computa- 
tional effort is to determine the time of the next collision 
with the oscillating bar, where the smallest root of an 
implicit equation (with possibly many roots) has to be 
found, see e.g. Refs. [30|, |3l|. The main quantity of inter- 
est is the time-evolution of the ensemble averaged mod- 



ulus of the velocity in y-direction (v x 
which is given by 
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where Vi(t) is velocity of the i-th particle at time t. The 
results of the simulations are shown in Fig. [2] for 3000 
oscillations of the bar on a semilogarithmic scale. The 
ensemble averaged velocity clearly grows exponentially 
(H)(£) ~ exp(Rt) as reported in Ref. 29], with a growth 
rate of R « 1.1 x 10 -5 . Let us now develop a physical 
picture of this acceleration process, answering the ques- 
tion how we can link the microscopic dynamics of single 
trajectories to the appearance of exponential Fermi ac- 
celeration. 



III. CONNECTION WITH THE FERMI-ULAM 
MODEL 

Since the v x component of a particle's velocity stays 
constant, see Fig. [TJ as long as the particle is in one of 
the interaction areas, the dynamics in y direction cor- 
responds exactly to the one of a particle moving in an 
one-dimensional (ID) Fermi-Ulam model (FUM), where 
the distance between the equilibrium position of the mov- 
ing wall and the static wall is given by h = H/2. The 
time ti the particle spends inside the interaction area is 
simply given by tj = l/v x . We define the time i^ the 
particle spends in the FUM by each passing of the in- 
teraction area (A or B, see Fig.Q} as the time difference 
between the first and the last collision with the oscillat- 
ing bar while the particle is in the interaction area. The 
motivation for this definition is that only collisions with 
the oscillating bar change v v , i.e. we want to keep track 
of the time in which a certain change in the velocity takes 
place. The times ti and £f are not identical, since once 
a particle enters the interaction area, a certain amount 
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FIG. 2. (Color online) Semi-logarithmic plot of the time- 
evolution of the ensemble averaged modulus of the velocity 
(|«|)(i). The velocity grows exponentially, {\v\)(t) ~ e Rt , with 
a growth rate of R = 1.1 x 10 5 (a = 0.1, w = 0.02). 




FIG. 3. (Color online) Phase space of the one-dimensional 
Fermi-Ulam model. For low v there is a chaotic sea, whereas 
for high velocities there are invariant spanning curves. The 
difference between the maximum and the minimum of such a 
curve grows with increasing v. The red lines show the analytic 
results for the invariant spanning curves based on the static 
wall approximation (SWA). In the inset, a single invariant 
spanning curve is shown for different driving amplitudes a. 
For large a, the SWA deviates significantly from the exact 
result. 



of time will elapse, before it collides with the bar. How- 
ever, for high velocities v, which means that v y is large 
since v x = const., tp converges towards tj. By 'high' we 
mean that v is large compared to the maximal velocity 
of the bar, i.e. for \v\ 3> wao (which implies \v\ ~ \v y \) 
we obtain tp « ti = l/v x . 



Since the dynamics of particles can be described for 
some time spans tp as a ID FUM, it is convenient to 
summarize some of the properties of the FUM (for a more 
detailed description see Ref. [18| and references therein) . 
The phase space of the FUM is shown in Fig. [3] For low 
velocities there is a large chaotic sea containing many 
regular islands (this regime is shrinked to a narrow band 
< v < 0.3 in Fig. \S§. Above the first invariant span- 
ning curve (FISC) with velocity v c , the motion becomes 
more and more regular, until for v 3> v c there are ex- 
clusively invariant spanning curves corresponding to a 
synchronized motion between the oscillating wall and the 
particles. Due to these invariant curves, there is no diffu- 
sion in momentum space and the FUM with a harmonic 
oscillation of the wall does not show Fermi acceleration. 

The invariant spanning curves v lsc ((f>) are not just 
straight lines, but show a characteristic shape (see 
Fig. [3]). There are infinitely many of them, which can be 
labeled by the velocity v lsc (cf> = 0) and are parametrized 
as v lsc = v lsc (cf>, v). The minimum of these curves is al- 
ways at <j) — and the maximum at (f> — tt/uj. The dif- 
ference Av lsc 



, ,1SC / 



v"~(0,v) grows linearly with 
v. For high velocities (v ^$> v c ) 



v lsc (n,v 
increasing v, i.e. Av mc «~ 
this can be rigorously shown within the so-called static 
wall approximation (SWA) [Jj| [l9j , which assumes that 
the oscillating bar is fixed in coordinate space but trans- 
fers momentum as if it would be moving. The distance 
between two collisions is then simply 2h and the time be- 
tween two collisions is At = 2h/v, where v is the velocity 
after the preceding wall collision. After a collision with 
the (only in momentum space moving) wall, the velocity 
Vi of the particle is 



v 1 = v isc (t,v)-2v w (t), 



(2) 



where v w (t) = yb(t) — — awsin(u;£) is the velocity of the 
wall. Since for high velocity the particle moves on an 
invariant curve we set v\ = v lsc (t + At, ii) and get 



v isc {t + At, v) = v isc {t, v) - 2v w (t). 



(3) 



Expanding this equation into a Taylor series up to first 
order and applying the continuous limit At — > 0, we ob- 
tain 2hv lsc (t,v)/v lsc (t,v) — — 2v w (t). Integration yields 

w isc (^,{i) = {)e (o/h) ' [1 - cos ^ 0g[O,2tt) (4) 

Therein, <j) is the phase of the wall oscillation, <j) = cut 
mod 2-7T. Obviously, the difference Av lsc between the ve- 
locity maximum and minimum of an invariant spanning 
curve is 



Av lsc = 



v isc {n, v) - v isc (0, v)=v- (e 2a/h - l) 



(5) 



and thus proportional to v. The v lsc (4>, v) of Eq. ([JJ 
are shown in Fig. [3] as red lines. We obtain a very good 
agreement with the exact results from the numerical sim- 
ulations. However, the inset shows a single invariant 
spanning curve for different driving amplitudes a, and 
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FIG. 4. (Color online) A typical trajectory in the investigated 
setup (see Fig. [1]). The particle enters four times the inter- 
action area (curves 1,3,4 area B, curve 2 area A). For a fixed 
time £f (or phase A(/> = uitF mod 2ir) it moves along the 
invariant spanning curves of the corresponding Fermi-Ulam 
model (thin red lines) , see also Fig. [3] It then leaves the in- 
teraction area with a certain velocity before it re-enters area 
A or B with the same velocity but now with a different phase. 



for large a, the SWA (Eq. 01 deviates significantly from 
the exact result, i.e. the SWA is valid for small driving 
amplitudes only. 

Based on the above phase space description of the 
FUM we can investigate the microscopic dynamics of a 
typical trajectory moving inside the oscillating bar bil- 
liard. Since we are ultimately interested in the acceler- 
ation process, we can assume that the particles are fast 
compared to the motion of the bar, i.e. again v ^> aau. 

This means that the velocity is also much larger than 
the velocity v c of the FISC. A particle enters the interac- 
tion area, let's say above the bar (interaction area A), at 
time t\ with a certain velocity v\ = \v\\ and spends the 
time tj « tp in it. Since it is fast, it collides many times 
during the time tp with the bar, thus moving along the 
corresponding invariant spanning curve of the FUM and 
leaving the interaction area at time t 2 = ti + tF with ve- 
locity vi . Now the particle propagates in the free part of 
the rectangular billiard (i.e. the part where no collisions 
with the bar take place) for a time % — 2d/v x . During 
this time, the modulus of the velocity does not change 
and thus the particle re-enters at time £3 = £2 + tb the 
interaction area with velocity v 2 . The particle can enter 
either the interaction area A or B, i.e. above or below 
the bar depending on the exact dynamics. However, for 
high velocities, in a good approximation this is a random 
process, as argued in Ref [29|], i.e. the particle will be 
injected with probability one half above and with prob- 
ability one half below the bar. If the particle is injected 
in part A, it re-enters the same FUM as described above, 



now at the phase 4> = ujt^ mod 2ir. If it is injected in 
part B we have to add a phase shift of 7r, i.e. <j) = (ut^+Tr) 
mod 2-7T, since the two FUM's (above and below the bar) 
can be transformed into each other simply by shifting 
the phase by it. Again, the particle moves for a time tp 
along an invariant curve and the whole described pro- 
cess - random injection in the FUM, leaving the FUM, 
re-entering etc. - starts over again. Exemplary, this pro- 
cess is shown for a typical trajectory in Fig. |4j The thick 
black lines show the motion along the invariant spanning 
curves and thus the evolution of the velocity of the parti- 
cle for four such injections into the interaction area. The 
thin red lines show the invariant spanning curves of the 
corresponding FUM's. We see perfect agreement, sup- 
porting our assumption that the particle moves for some 
time tF along such invariant curves. Note that in the 
above terminology, the curve 2 corresponds to the FUM 
of the interaction area A, whereas the curves 1,3,4 corre- 
spond to the phase shifted FUM of the interaction area 
B. 

Let us describe the above process in a more quantita- 
tive way. The particle enters the interaction area at t\ 
with a high velocity v\ 3> a^ui and v\ 3> v x . The corre- 
sponding invariant spanning curve on which the particle 
will move for the time tF can be calculated as follows: 
The entry phase (j}\ is given by 4>i = uit\ mod 27r; by 
setting v x = v isc (c/>i,v) = {^AH 1 - 008 ^] we obtain 
v = Ul e(- Q o/' l )-[ 1 - cos< # , i]. The exit velocity is then (re- 
member tp = l/v x ) 



v 2 = v ■ 



+ ul/v x ,v) = V\ 



gf[l— cos(0i +uil/vx)] 
e f[l-cos0i] 



(0) 



This procedure can be repeated again and again, yield- 
ing v n = X)j=i ^i + w ii where the Avi are obtained 
in the above described manner by exploiting the piece- 
wise motion on the invariant spanning curves. However, 
this sum is not suitable to obtain a closed expression for 
v(t) (or (H)(£)), especially since it contains the random 
phase shifts of it. We thus employ in the next section a 
random walk model, based on the statistical properties 
of the just described procedure, and will thus be able to 
calculate explicitly the exponential growth rate. 



IV. RANDOM WALK MODEL 

The velocity in the random walk model is written as: 

v n+ i = v„ + Av n , (7) 

where Av n can be positive or negative. Since the Au„ are 
determined by moving along parts of the invariant span- 
ning curves of the FUM, and for the latter we know that 
Av isc ~ v, we conclude that the Av n are proportional to 
v n , so Av n = ±cv n , yielding 



Vn+i = v n ±c- v n = (1 ± c)v n 



(8) 



The constant c is an effective constant, besides the ge- 
ometry of the billiard, c depends in particular on the 



entry phase <f> of a particle into the FUM and of course 
on the time tp — l/v x (or phase A<j> = ujI/v x mod 2-k) 
the particle spends in the FUM. We consider an ensem- 
ble of particles all starting with the same v x , thus the 
phase shift Acf> is the same for all particles and A(j> is a 
constant. Under these assumptions we will proceed as fol- 
lows: • We show that the random walk model of Eq. ((5J) 
leads to exponential Fermi acceleration. • We will de- 
termine the effective c using the phase space properties 
of the FUM. To this end, we consider a particle after N 
steps, i.e. N cycles through the FUM. The probability 
p(k) to have completed k more positive than negative 
steps Aii is then 



p(k) 



1 



N 



2 N \(N + k)/2 



(9) 



We define v^k) as the velocity that is reached after N steps 
with k steps more in positive than in negative direction 
and U(-fc) as the velocity that is reached after TV steps 
with k steps more in negative than in positive direction. 
There are of course many different paths, leading to the 
same vha, however the order of the steps is irrelevant. 
Let us assume (without loss of generality) N to be even, 
then k can be any even number between —N and N. The 
expectation value of the modulus of the velocity after 
N steps is given by summing over all possible A;'s and 
weighting them with the corresponding probability p(k): 



<M> 



N 

N= Yl P( k )\ l '\(k) 

k— — N, k even 

N/2 



= - y 



2 N ^ V JV/2 + k 



N 



M(2fc) 



(10) 



By setting 7 := 1 + c and vo as the initial velocity, using 
Eq. © for the positive sign, we obtain 



<M> 



N 



Vq 



(1 + 7 2 ) 
(27)" 



N 



(11) 



To switch from the number of cycles N to the actual time 
t, we use that the time between two collisions with the 
vertical walls is given by L/v x . Substituting N = v x t/L 
yields the exponential time-law 



Rt 



(\v\)(t) = v e 
where the growth rate R is given by 



R= JL m ( 7 / 2 + I/27). 



(12) 



(13) 



The random walk model reproduces an exponential de- 
pendence of the ensemble averaged velocity on time. 
Nevertheless, in order to determine the system-specific 
value of the growth rate R, we still have to determine 
the effective constant c. To this end, we rewrite the part 
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FIG. 5. (Color online) Frequency dependence of the growth 
rate 7?. The red line shows the result from the numerical 
simulations and the blue line the growth rate as predicted by 
our random walk model, c.f. Eqs. (|13|) and (|17[) (a — 0.025). 
The inset shows the amplitude dependence (w = 0.1) of the 
growth rate. 



of Eq. IjHJ) with the '+' as c = (v„+i — v n )/v n . The ve- 
locity i> n +i is, according to Eq. ©, given by 



V n +1 



+ Atp,v(v n ,<pn)), 



(14) 



where we write v = v(v„,(j> n ), since v depends on 
the entry phase </> n and the entry velocity v n . For a 
fixed A(f), we thus have c = c(v n ,(f> n ) — v lsc (v n ,4> n + 



1 with v" 



oc v n 



leading to 



A(j))lv lh -[v n ,(p n (, 
c{v n ,4> n ) — c((f> n ). Thus, the effective c is given by av- 
eraging over all entry phases <p n that lead to a positive 
Av n (This is sufficient, since the random walk model of 
Eq. (fg|) , (fTT)]) has intrinsically included the '-' part, allow- 
ing the Av to be negative.): 



Ceff 



l>n,2 



C(tj> n )d(j) n , 



where 



c(4> n ) 



Qh 



[l-cos(<A»+A0)] 



1. 



(15) 



(16) 



The integral over <p n has to be evaluated such that Av n is 
positive in order to account all accelerating trajectories, 
i.e. 4>n,i = —Acf)/(2uj) and § n ^ — ir/u — A0/(2w) yield- 
ing for the normalization A^ = l/(<f) n ,2 — <t>n,i) — w/vr. 
Since the growth rate R depends on 7, see Eq. (fl~3|) . and 
7 = 1 + c, we finally obtain 



7eff = 1 + C e ff 



7r/w-A0/2w 
LU f gf [l-cos(0„+A0)] 



-A<t>/2u 



1 



<C), 



(17) 



Here we assumed that the entry phases <p n are uniformly 
distributed, which is valid on a sufficient long time scale if 



the ratio L/v x is incommensurate with the driving period 
T = 2tt/uj. Note that the growth rate R of Eq. JT5J to- 
gether with the corresponding 7 of Eq. (flTf has been ob- 
tained ab initio without any fit parameters. However, the 
corresponding values are too small. From the results of 
the simulation shown in Fig. [2j for a = 0.1 and lj = 0.02 
we obtain for the growth rate R f=a 1.1 ■ 10~ 5 . Inserting 
Eq. HID into Eq. ([13]) leads to R = 0.51 • 10~ 5 , which 
provides the order of magnitude of the numerical result, 
but is too small by about a factor 2. The reason for this 
is as follows: The sum of two sequences of the dynamics 
each consisting of e.g. three steps, one with a c((j)) close 
to the maximally possible value c max in each step and 
one with steps close to the minimal value of c(</>) (i.e. 
Cmin = 0) contributes more significantly to the ensem- 
ble average (\v\), than the sum of two corresponding se- 
quences both with step width c = (c max + c m ; n )/2. How- 
ever, according to Eq. (JT7J) we calculated a mean of the 
latter type and therefore obtained a lower bound of the 
correct growth rate. The deviation of the result obtained 
from the simulation and the one from the model can thus 
be understood as a consequence of the negligence of cor- 
relations. One way of effectively including correlations 
to the definition of c c g is to define c c g = (c m ) / '(c™ 1-1 ) 
with m > 1. Therein (c m ) , denotes the average of c m (0) 
over all phases in the interval [— A0/(2w), (tt — A0/2)/w] 
with c given by Eq. (TT6|) . An upper bound for the growth 
rate can be obtained within the assumption that all steps 
but those with maximal c(0) are suppressed, i.e. by cal- 
culating the effective c with the assumption that c(</>) is 
equal to the maximally possible value c max at each step: 

e f[l-cos(7r/2+A0/2)] 
7max — -L t C m; 






[1-cos(tt/2-A0/2)] 



1 



lim 

N—xx 






(18) 



This upper bound leads to R = 1.24-10 -5 , which is quite 
close to the result obtained from the simulation. In order 
to test these estimations for a whole range of parameters 
and also to show that our random walk model correctly 
describes the whole dependency of the growth rate of the 
parameters of the system, we extract the growth rate 
R for different driving frequencies lj at fixed amplitude 
a = 0.025 by performing a numerical simulation for each 
value of the frequency and compare R with the corre- 
sponding result obtained from our random walk model. 
The growth rate R{uj) (see Fig. [5]) shows characteristic 
(decaying) oscillations, as already theoretically predicted 
in Ref. [29(. The minima where R(lj) is exactly zero can 
be easily understood. At these values of ui, the driving 
period T = Itx /lj and the time between two collisions 
with the same vertical wall 2L/v x are commensurable, 
what leads to a 0-periodic entering and leaving of the 
FUM for v y — > 00 (when the first and the last collision 
with the oscillating bar are converged to its edges) . Note 
that the occurrence of the minima in the growth rate 
R(u>) are based on the fact, that all particles of the en- 
semble possess the same, constant velocity in x-direction. 



Apparently, these characteristic oscillations are fully re- 
produced by our model. 

The inset of Fig.[S]shows the analogous comparison be- 
tween model and simulation for a fixed driving frequency 
lj = 0.1 and different but small values of the amplitude 
(the regime where the invariant spanning curves of the 
FUM can be well approximated within the SWA). From 
this, we firstly observe that the growth rate strongly in- 
creases with the amplitude of the oscillating bar, and 
secondly, that this dependency can be well explained by 
the random walk model. According to the good agree- 
ment between the simulation and the model, we may 
conclude from Eq. (fTTj) and Eq. (fT3[) that the amplitude- 
dependence of the growth rate is approximately given 
by R(a) oc ln(cosh(a//i)). As Fig. [3] reveals, we may 
not expect that this is also true for large values of the 
amplitude, since then the expressions for the invariant 
spanning curves obtained within the static wall approxi- 
mation strongly differ from the numerical results. Obvi- 
ously, the result of the simulation is between the estima- 
tions for the lower and the upper bound for all values of 
the system parameters. 

These results indicate that all details of the specific sys- 
tem under consideration which are not accounted for the 
random walk model including the existence of a chaotic 
sea do not contribute crucially to the growth rate. Even 
more, all details of the specific system are only needed 
to calculate the effective c. The requirements for the oc- 
currence of exponential Fermi acceleration are compara- 
tively weak: A temporally periodic entering and leaving 
of invariant spanning curves with t^ ax 
different phases is sufficient. 



V. CONCLUSION 

In this work, we have investigated the physical mech- 
anism leading to exponential Fermi acceleration (FA) in 
the rectangular billiard with an oscillating bar inside. In 
particular, we showed that the dynamics of individual 
trajectories can be understood as alternating phases of 
motion in an appropriate one-dimensional Fermi-Ulam 
model (FUM) and free propagation. During the tempo- 
ral FUM phases, the particles move (in the high velocity 
regime which is of interest here) on invariant spanning 
curves of the FUM, which can be - at least for small driv- 
ing amplitudes - obtained analytically within the static 
wall approximation. Using the intrinsic property of the 
invariant spanning curves of the FUM that the difference 
between the maximal and the minimal velocity grows lin- 
early with the minimal velocity of the invariant curve 
and the fact that acceleration and deceleration have equal 
probability, the process can be modeled as a random walk 
with step width proportional to the velocity itself, i.e. 
w n +i = v n ± cv n . This model explains the occurrence 
of the exponential acceleration. Calculating an effective 
step width e e g we obtain a good ab initio estimation of 
the growth rate and reproduce the whole qualitative de- 



pendency of the system parameters. We emphasize, that 
our random walk model reflects that a temporally pe- 
riodic entering and leaving of equally shaped invariant 
spanning curves which have the property that the dif- 
ference between the maximal and the minimal velocity 
grows linearly with the minimal velocity of the latter 
are the key ingredients for the occurrence of exponential 
Fermi acceleration (all details of the specific system are 
contained in the factor c). This opens the perspective of 
searching for other systems exhibiting the phenomenon 



of exponential Fermi acceleration. 
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